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Abstract 



1-^ ■ We analyse hysteresis in a one-dimensional anti-ferromagnetic ran- 

Ph ', dom field Ising model at zero-temperature. The random field is taken 

to have a rectangular distribution of width 2A centered about the 
origin. A uniform applied field is varied slowly from — oo to +oo and 
back. Analytic expression for the hysteresis loop is obtained in the 
^ ' case A < \J\, where | J| is the strength of the nearest neighbor inter- 

j^ ■ action. 
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1 Introduction 



In a recent paper ^ , we presented an analytic solution of 
the zero-temperature dynamics of the anti-ferromagnetic ran- 
dom field Ising model in a slowly varying uniform applied field 
over a truncated range of the field. The purpose of the present 
paper is to extend the earlier solution to cover the entire range 
of the applied field, and hence to obtain an analytic expression 
for the hysteresis loop. We refer the reader to reference [1] for 
a detailed discussion of the model, and the method of its solu- 
tion. Here, we recapitulate the salient points of reference [1] for 
maintaining continuity with it, and to set up the notation. 
The Hamiltonian of the system is given by 

H = -JJ2s^Si+i-Y.hiSi-haJ2sz, (1) 

i i 

where s^ = ±1 are Ising spins at sites i = 1,2,3... of a one- 
dimensional lattice, ha is a uniform applied field, and hi is a 
quenched random field drawn independently from a continuous 
bounded distribution, 



p{hi) =2^ if -A<h,<A 

= 0, otherwise. (2) 

The external field ha is slowly increased from — oo to +oo, 
and the system is allowed to relax by the single-spin-flip Glauber 
dynamics at zero temperature, i.e. spins are allowed to flip one 
at a time if it lowers their energy. This dynamics takes the 
system to a stable state where each spin has the same sign as 
the net local field /,• at its site. 



Si = sign /, = sign [J(s,_i + s,+i) + /i^ + ha] (3) 

The magnetization m{ha) per spin at the apphed field ha is 
given by 

miha) = -T.s^ (4) 

Our object is to calculate m{ha) for all values of ha in the 
range [—00 to +00]. In reference [1], we obtained m{ha) for 
[-00 < ha < 2|J| - A] in the case [A < \J\]. The result is 
expressed conveniently in terms of variables h and p defined as 
follows: 

h = -ha mod |2J|; p{h) = J^ p{h^)dh^ = 

The field h serves as a reduced applied field whose importance 
is explained below. If A < | J|, the spins turn up in three cate- 
gories. The first to turn up are the spins whose nearest neighbors 
are both down (ramp-I), then the spins with one neighbor up and 
one down (ramp-II), and finally the spins with both neighbors 
up (ramp-III). On each ramp, a spin with quenched field hi turns 
up when the reduced field h crosses the threshold h + hi = 0. 
Thus the reduced field acts as an effective applied field on each 
ramp. The reduced field covers the same range as the quenched 
field, i.e. [—A < {h,hi) < A]. The parameter p is the fraction 
of sites in the system with hi > h. This parameter also serves as 
a convenient measure of the reduced field on each ramp. A key 
quantity on each ramp is the fraction of spins whose quenched 
field is larger than the reduced field h. This quantity determines 
the shape of the ramp. The difficulty in calculating this quan- 
tity comes from the fact that the distribution of quenched fields 



on each ramp is modified a posteriori by the dynamics of the 
system. In reference [1], we obtained the following results: 



Ramp-I: [-2| J| - A < /i„ < -2| J| + A)] 



m{h) = —e 



-Mh) 



Plateau-I: [-2| J| + A < ha < -A] 

m{h) = e^ 

Ramp-II: [-A<ha<A] 

m{h) = - (3 + 8e"^ + e"^) + 8 (l + e"^) p - 4p'^ 



+ 



1 + e 



-1 



-1\ ^2 



8p-2[l + e-'}p' + -pHe-P 
- (5 + 2p) e-'^P 



Plateau-II: [A < ha < 2\J\ - A] 

■27 7 



m{h) = 



30 6 3 



Reference [1] did not contain results for ramp-III, i.e. m{ha) 
in the range [2| J| — A < ha < 2| J| + A]. This calculation is 
taken up in the following. 



2 The Ground State 

Plateau-II forms the ground state for the evolution of ramp-III in 
an increasing apphed field. The ground state has solitary down 
spins dispersed in a sea of up spins. Following the nomenclature 
of reference [1], we call the down spins singlets. The net field 
on a singlet is equal to —2\J\ + hi + ha- This lies in the range 
[-2\J\-A+ha] to [-2|J| + A + /iJ. Therefore, the singlets begin 
to turn up at ha = 2\J\ — A, and are all up at ha = 2| J| + A. 
These limits mark the beginning and the end of ramp-III. The 
shape of ramp-III is determined by the fraction of up spins at 
applied fields between these limits. We take an arbitrary value 
of the applied field ha = —2\J\ — h on ramp-III, and calculate 
the fraction of singlets which are up at this field. In other words, 
we calculate the fraction of singlets with hi > h ot equivalently 
Pi < j9, where 

A-h, 
Pi = 



2A 

As discussed in reference [1], the singlets on plateau-I fall 
into three categories. In category (1) are the singlets which 
were created on ramp-I. The fraction of singlets which belong 
to this category, and whose quenched field is greater than h is 
given by, 

Pi'ip) =P-h^- e-''] - - [e-' - (1 - P)] (5) 



2 



e 



In category (2) are the singlets which were created on ramp-II 
by a vanishing doublet. The fraction of singlets in this category 
whose quenched field is greater than h is equal to. 



P^\p) = [e-P-{l-p)Y (6) 

Finally, category (3) singlets are those that were created on 
ramp-II by the central spin of an unstable up triplet flipping 
down. The fraction of singlets in this category whose quenched 
field is larger than h is given by, 



-1 



l + 2e-M + 



Pi'ip) 



- [l + 5e-^ + 4e"^] p^ + - 



1 



-1 



+ 2e"' + 4e 



-2 



- + 4e-i + e-2 



p 



1 

4 



1 
20^ 



1 + er^\ / + -^/ + [l + 2e"^ 



.2^e^(i+p)+_p2e-^+- 



1 



^2p 



(7) 



The total fraction of singlets present on plateau-II whose 
quenched fields is larger than h is equal to. 



P'\V) = P('{P) + Pf{p) + Pf{p) (8) 

The above equation gives the fraction of singlets which will 
turn up in an applied field ha = 2\J\ — h from among the singlets 
initially present on plateau-II. However, when these singlets turn 
up, some of their nearest neighbors turn down. This process 
creates new singlets [0]. The newly created singlets turn up 
later at a higher applied field on ramp-III. We need to calculate 
the fraction of newly created singlets, and their restoration to 
the original state as function of the applied field before we can 
determine ramp-III. 



3 Creation of New Singlets 

In this section, we consider the circumstances in which the de- 
struction of a singlet on ramp-III accompanies the creation of a 
new singlet. Consider a singlet in the ground state, say the spin 
at site 3 in Figure 1. Its nearest neighbors at sites 2 and 4 are 
up. A next nearest neighbor at site 1 is down, and the other 
next nearest neighbour at site 5 is up. Suppose the singlet just 
flips up on ramp-III at an applied field ha, i.e. hs — 2\J\ + ha = e, 
where e > 0. The applied field at this point is /ia = 2| J| — /i3 + e. 
We ask the question, could a nearest neighbor of the singlet flip 
down when the singlet flips up? 

First, consider the nearest neighbor at site 2 in Figure 1. 
After site 3 has flipped up, site 2 has one nearest neighbor up 
and one down. Thus site 2 may flip down if h2 + ha < 0, or 
h2 ^ hs — 2\J\ — e. However, this is not possible because /i2 and 
/i3 lie in the range [—A to +A], and we are considering the case 
A < \J\. Thus a nearest neighbor of a singlet which has both 
its nearest neighbors down will stay up when the singlet turns 
up as long as A < |J|. 

Next, consider the spin at site 4. After the singlet has turned 
up, site 4 has both its nearest neighbors up. It will turn down 
if -2| J| + /i4 + /ia < 0, or hi < hs - e. If h^ < h^, site 4 will 
turn down when site 3 turns up. After site 4 has turned down, 
site 3 has one neighbor up and one down. The net field at site 
3 is then /i3 + /z-a = 2| J| + e, which is positive. Thus site 3 will 
stay up after site 4 has turned down. We conclude that, when 
a singlet turns up, its nearest neighbor may turn down if that 
nearest neighbor has less quenched field than the singlet, and 
also if it had one nearest neighbor already up before the singlet 



turned up. Note that when site 4 turns down, it increases the 
upward field at site 5. There is no possibihty of site 5 turning 
down as a result of site 4 turning down. Consequently, a spin 
turning up on ramp-III does not cause any change in the state 
of spins beyond the nearest neighbors (absence of avalanches). 
Before proceeding further, we rewrite the two rules which will 
guide us in the following analysis. 

Rule-1: When a singlet turns up on ramp-III, its nearest neighbor 
stays up if the adjacent next nearest neighbor is down, and 

A<|J|. 

Rule-2: When a singlet turns up on ramp-III, its nearest neighbor 
turns down if both of the following conditions are satisfied: 

(a.) the adjacent next nearest neighbor is up, and 

(b.) the nearest neighbor has less quenched field than the 
singlet. 

The next question is, in what circumstances, the quenched 
field on a nearest neighbor of a singlet on plateau-II can be 
smaller than the quenched field on the singlet. We examine var- 
ious possible cases. Suppose the singlet in question was created 
on ramp-II by the destruction of a doublet. For example, look 
at Figure 1 again and suppose sites 2 and 5 were up, and 3 and 4 
were down on plateau-I. This requires h^ < h2 and h^ < h^. Let 
hs < h^ without any loss of generality. In this case, site 4 will 
turn up on ramp-II. When site 3 turns up on ramp-III, no new 
singlet can be created because both nearest neighbors of site 3 
have a larger quenched field than site 3 (Rule-2). 

Next, let us assume that the singlet in question was created 
on ramp-II by the unstable central spin of an up triplet fiipping 



down. The central spin of an up triplet flips up for the the first 
time on ramp-I. It separates two adjacent doublets on plateau- 
I. Therefore, the quenched field on the central spin is larger 
than the quenched field on each of its nearest neighbors. The 
nearest neighbors of the central spin flip up on ramp-II (the 
central spin flips down at this event). Therefore each nearest 
neighbor of the central spin has a larger quenched field than the 
next nearest neighbor of the central spin which is adjacent to it. 
Consequently the central spin has a larger quenched field than 
each of its next nearest neighbors. When the central spin flips up 
for the second time on ramp-III, its next nearest neighbors are 
still down. Therefore there is no chance (Rule-2) for the nearest 
neighbors of the central spin to flip down when the central spin 
flips up on ramp-III. 

Having ruled that the destruction of a singlet created on 
ramp-II does not give rise to a new singlet on ramp-III, we are 
left with singlets created on ramp-I and present on plateau-I. 
When one of these singlets disappears on ramp-III, it may cre- 
ate a new singlet if one of its nearest neighbors has a smaller 
quenched field than the singlet, and if the next nearest neighbor 
is up. We note that only one of the two nearest neighbors of a 
singlet on plateau-I can have a quenched field which is smaller 
than the singlet. A site which has a larger quenched field than 
both of its nearest neighbors must be up on plateau-I. Now con- 
sider a singlet on plateau-I as shown at site 2 in Figure 2. Let its 
nearest neighbor on the right (at site 3) have a smaller quenched 
field (/i3 < /i2). Site 4 must be down because there are no strings 
of up spins of length greater than unity on plateau-I. Site 5 may 
be up or down. The two possibilities for site 5 are depicted in 
Figures 2 and 3 respectively. In Figure 2, a singlet is followed 



by a singlet on plateau-I. In Figure 3, a singlet is followed by a 
doublet. 

Consider Figure 2 first. What is the probability per site that 
such an object occurs on plateau-I? A singlet on plateau-I must 
be followed by a singlet or a doublet. It was shown in reference 
[1], that the probability per site of the occurrence of a doublet 
is equal to e^^, and the probability that a doublet is followed 
by a doublet is equal to |e^^. Therefore, the probability that a 
doublet is followed by a singlet (or vice-versa) is equal to |e"^. 
It was also shown that the probability of the occurrence of a 



singlet is equal to ^ 1 — 3e 



Therefore, the probability per 



site that a singlet is followed by a singlet on plateau-I is equal 



to i [l - ¥e 



2 



13^-2 



3 



We have obtained above the probability of the occurrence 
of objects shown in Figure 2. Our immediate interest lies in 
a subset of these objects with h^ < min (/i2,/i4). This subset 
determines the creation of new singlets. Suppose /i2 < h^. Then 
on ramp-III, site 4 will flip up before site 2. When site 2 flips 
up, site 3 will flip down because the conditions for the creation 
of a new singlet are fulfilled (Rule-2). In order to calculate the 
fraction of newly created singlets, we need to know the distri- 
bution of fields /i2, /is, and h^. These may be obtained on the 
lines of reference [1] if one notes that sites 1 and 5 must have 
flipped up on ramp-I before site 3. Site 1 must have flipped up 
before site 3 because hi > h2 > h^. Similarly, site 5 must have 
flipped up before site 3 because h^ > h^, h^ > h2, h2 > h^, and 
therefore h^ > h^. This proves that just before site 3 flipped up 
on ramp-I, sites 2, 3, and 4 formed a string of three down spins 
bordered by up spins at sites 1 and 5. The screening property 
discussed in reference [1] can be applied here to conclude that 

10 



the distributions of /i2, and h^ are independent of each other, 
and each is given by 

Prob[p < P2 ^ P + dp] = [1 — e^^]dp 

Prob[p < P4 < p + dp] = [1 — e^^]dp 

where 

A-h, 
Pi = 



2A 

Given that /i2 hes in the range [h to h+dh], h^ must be uni- 
formly distributed in the range [—A to h]. Thus, 

Prob[p < ps < p + dp] = dp (if p^ > P2) 

The contribution of Figure 2 to the fraction of newly created 
singlets when all sites with hi > h have been relaxed on ramp-III 
is given by, 

Pl''{p) = 2 1' [1 - e-P^] dp2 £ dps I' [1 - e~P'] dp, 

= 5 _ 2p{l-p) - ^p' - 2{l-p + p')e-P+^{l - 2p)e-'P 

(9) 

Consider Figure 3 next. It shows a singlet at site 2 followed 
by a doublet at sites 4 and 5. An identical contribution will 
come from a situation where the doublet occurs at sites 2 and 3, 
and the singlet at site 5 as shown in Figure 4. We work out the 
contribution from Figure 3 explicitly, and then multiply it by 
a factor 2 to take into account the contribution from Figure 4. 
When site 2 flips up, site 3 will flip down if h^ < /i2, and if site 
4 is up as well. Note that site 4 would have flipped up earlier on 

11 



ramp-II if h^ > h^. If /is > /i4, site 5 would have flipped up on 
ramp-II. In this case there is no possibility of site 3 flipping down 
when site 2 flips up. The reason is as follows. The presence of 
a doublet at sites 4 and 5 on plateau-I means that h^ < h^ and 
h^ < Hq. We must have h^ < h2 for the creation of a new singlet 
(Rule-2). If /i4 < /i3, and h^ < /i2, then we have h^ < h2. Thus, 
when site 2 turns up on ramp-III, site 4 will be down, and one 
of the two conditions for the creation of a new singlet is not 
satisfied. 

A similar reasoning as applied in the analysis of Figure 2 
reveals that site 3 must have flipped up after sites 1 and 6. The 
distribution of fields at sites 2 and 5 in Figure 3 must be similar 
to the distribution of fields at sites 2 and 4 in Figure 2. The 
distibution of fields at sites 3 and 4 must be uniform over the 
interval [—A to min(/ii, ho)]. Thus the contributions of Figures 
3 and 4 to the fraction of newly created singlets when sites with 
hi > h have been relaxed on ramp-III is given by, 

Pf'ip) = 2 r [1 - e-P^] dp2 /' dps /' dp, /' [1 - e-P'] dp5 
-^/ + { (^ + 3e-i) - (1 + 2e-i) p + e" V - ^-p'^ 



-o" > e P- e-^P 



(10) 



4 Destruction of New Singlets 

The destruction of newly created singlets on ramp-III can be 
analysed in a similar manner as their creation. For example, 

12 



refer to Figure 2 again, and recall that /i2 = min(/i2, /i4), and 
^3 ^ ^2- In this figure, the creation of new singlets when all 
sites with hi > h have been relaxed is determined by /i2 > h. 
The destruction of new singlets when all sites with hi> h have 
been relaxed is given by h^ > h. The result is an integral similar 
to the expression for the creation of singlets. Only the limits of 
the integration are altered. We obtain. 



)//// 



P3 r 



1 

2 



-P2 



^p.i:ii 



-PA 



1 



-2p 



dp4 

2pe'P + p{l -p) + -p^ 

o 



:ii) 



Similarly, the contribution from Figures 3 and 4 is given by. 



pr{p) = 2dpsi 1 



P3 



-P2 



dp2 / dp4 / 1 — e ^^ dpr, 

Jpz Jpi '- ■' 



-M/ 



- { (1 + 2e-i) -2(1 + e-i) i? + /} e"^ + e-^P 



1 + 36"') p'' + -/ 



(12) 



5 Ramp-III 

Putting the various terms together, the probability that a ran- 
domly chosen spin on the chain is up on ramp-III is given by 



Pi" 



(p) = p,"(i) + p("{p) + pi"(p) + Pi" 



(p) 



-Pl"(v) - Pi"{p) + Pi"(p) + Pi"(p) 



(13) 
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The magnetization on ramp-III is given by, 

m'''{h) = 2Pl''{p) - 1 (14) 

After some simplification, we obtain 

- (l - be-^ + 4e"2) / + - (l + 2e"^ + 2e-^) / 

- I (^ + lOe-i) - (2 + 4e-i) i> - (4 - 2e-') / - ^/} e"^ 

+ (4 + 2j9))e-2?(l5) 

The above expression has been superposed on the simulation 
data in Figure 5. The agreement is excellent. The analytic result 
is indinguishable from the simulation on the scale of Figure (5). 

6 Hysteresis Loop 

So far, we have analysed the magnetization m{ha) in increasing 
applied field. We have shown that the analytic result agrees 
with the simulation rather well. The magnetization mji{ha) in 
decreasing field (return loop) is related to m{ha) by a symmetry 
of the model, i.e. mR{ha) = —m{—ha). Thus we have implicitly 
determined the return loop as well. The return loop has been 
shown in Figure (5) by a broken line. The hysteresis in the 
anti-ferromagnetic RFIM is rather small, and the two halves 
of the hysteresis loop lie very close to each other. In order 

14 



to show them more clearly, we have plotted in Figure (6) the 
separation between the two halves of the hysteresis loop versus 
the applied field. To be precise, we have plotted [m{ha) — m] 
and [mR{ha) — m] vs ha, where m = ^ [Tnj{{ha) + 'm{ha)]- 

As we may expect, the agreement between the theoretical ex- 
pression and the simulation is excellent on the scale of Figure 

6 as well. In Figure 6 the simulation data was obtained from 
a system of 10^ spins, and averaged over 10^ independent re- 
alizations of the random field distribution. The set of applied 
fields where the spins fiip on each half of the hysteresis loop is of 
course different for each realization of the random field distribu- 
tion. Therefore, the average over different realizations requires 
a judgement on how to group the data. We divided the entire 
range of the applied field from [— 2|J| — A to 2|J| + A] into 
3 X 10^ sections (bins) of equal width. The data in each bin was 
averaged separately. The simulation data shown in Figure 6 is 
a much sparser set of data (in order not to crowd the figure). 
We have shown the simulation data at intervals of 6 ha = -1, and 
the theoretical expression at intervals of 6ha = .01 (joined by a 
continuous line on the lower half, and a broken line on the upper 
half of the hysteresis loop). 

7 Discussion 

We have considered the zero-temperature dynamics of a one- 
dimensional anti-ferromagnetic random field Ising model, and 
obtained an analytic solution of the model if the following con- 
ditions apply: 

1. All spins are down initially, and the applied field is swept 

15 



from ha = —oo to ha = +00 infinitely slowly. The solution 
is also applicable by symmetry to the case when all spins are 
up in the initial state, and ha is decreased from ha = +00 

to ha = — CXD. 

2. The random field has a uniform bounded distribution in 
the interval [— A to + A]. We considered the case A < \J\. 
The simplifying feature of this case is that the increasing 
applied field exhausts all strings of down spins of length 
three or more (ramp-I) before working on strings of down 
spins of length two (ramp-II). Similarly, strings of down 
spins of length one (singlets) are turned up (ramp-III) only 
after the doublets are finished. 

The second condition mentioned above has been adopted es- 
sentially for simplicity. It serves to illustrate the method of 
solution with a minimum of algebraic detail. We do not see a 
conceptual difficulty in applying the same method in the case 
A > I J|, or when the random field has an unbounded continuous 
distribution but we do not go into these details here. 

The restriction to an initial state where spins are either all 
down or all up appears to be necessary so far. A similar difficulty 
is encountered in the ferromagnetic random field Ising model 0] . 
The relaxation dynamics of the ferromagnetic model is qualita- 
tively different from that of the anti-ferromagnetic model. The 
relaxation process in the ferromagnetic model is abelian, while 
in the anti- ferromagnetic model it is not. We have made some 
progress in the methods of analytic solutions in both cases, but 
a qualitatively new idea appears to be needed in extending these 
methods to an arbitrary initial state. 
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Figure Captions 

Figure 1: 
A singlet (site 3) with one next nearest neighbor down (site 1), 
and one next nearest neighbor up (site 5). When the singlet 
turns up at an applied field ha, the spin at site 2 stays up if 
A < I J|, but the spin at site 4 flips down if /i4 < /is. 

Figure 2: 
Two adjacent singlets on Plateau-I: If /i2 = min(/i2, /i4), and 
hs < /i25 then the spin at site 3 will flip down when the spin at 
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site 2 flips up on ramp-III. This process creates a new singlet on 
ramp-III. 

Figure 3: 
A singlet followed by a doublet on Plateau-I: If h^ > h^, and 
^3 ^ ^2, then a new singlet will be created at site 3 when the 
spin at site 2 turns up on ramp-III. 

Figure 4: 
A doublet followed by a singlet on Plateau-I: If h^ > /i2, and 
hi < /i5, then a new singlet will be created at site 4 when the 
spin at site 5 turns up on ramp-III. 

Figure 5: 
Magnetization per spin m in an applied field ha for an anti- 
ferromagnetic RFIM (J = —1) for a rectangular distribution of 
the random field of width 2 A = 1. The analytic result in in- 
creasing applied field is shown by the solid line, and in decreasing 
field by a broken line. Numerical result from a simulation of 10^ 
spins averaged over 3 x 10^ different realizations of the random 
field distribution has been superposed on the theoretical result 
by dots. The numerical result is indistinguishable from the the- 
oretical result on the scale of this figure. 

Figure 6: 
Separation between the two halves of the hysteresis loop in Fig- 
ure 5 has been magnified by plotting it relative to the average 
value of the two halves. The solid line shows 6m = [m — m], 
and the broken line 6m = [ijir — rn] vs the applied field ha; 
m = I [niR + m] ; m and m.R are the magnetizations at applied 
field ha in increasing and decreasing field respectively. 




h 



3 



(D @ © © 



h, 



h 



4 



Figure 1 



^2 


i 


I 


K 


b ® 


© 


(4) (i 


' 




. . 


h 


■ 


' 



Figure 2 




h, 



h 



4 



h. 



h 



5 



@ ® © ® © 



Figure 3 




h, 



h 



3 



h 



4 



h 



5 



© © © © © 



Figure 4 



E - 



-2.5 



-1.5 



Figure 5 




-0.5 



h. 



0.5 



1.5 



2.5 



Figure 6 



0.02 - 



0.01 - 



f \ 



(P 



® 



O €> 0- €^0-0-0 -G -G <5 -€> -Q 



U 






.0" 



jO -€>-€>€>€> e^ €> 0^0-0-0- 



E 

OO 



-0.01 - 



-0.02 - 



d 




00OO0OO0OO 




OOOOOOOO0O 





h 



